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We compute and analyze couples of ground states of 3D spin glass systems 
with the same quenched noise but periodic and anti-periodic boundary condi- 
tions for different lattice sizes. We discuss the possible different behaviors of 
the system, we analyze the average link overlap, the probability distribution 
of window overlaps (among ground states computed with different bound- 
ary conditions) and the spatial overlap and link overlap correlation functions. 
We establish that the picture based on Replica Symmetry Breaking correctly 
describes the behavior of 3D Spin Glasses. 

I. INTRODUCTION 

Understanding the structure of the states of three dimensional (3-D) Ising spin glasses at 
finite temperature is a very interesting problem. Long time ago a mean field theory for spin 
glasses has been constructed and successfully compared with the numerical properties of a 
long range model, the Sherrington-Kirkpatrick model, where mean field theory is exact by 
definition [IJ. 

In the mean field approach one finds that the free energy landscape is strongly corrugated, 
and it is characterized by the presence of many local minima, which correspond to spin 
configurations very different one from the other. There are many low free energy local 
minima which have a total free energy which differs from the ground state free energy of 
an amount which is of order 1; these minima contribute to the probability distribution of 
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the spins also for arbitrary large volumes. In some sense one can say that in the infinite 
volume there are many equilibrium states (for a better qualification of the terminology see 
0): these states are locally different one from the other. For historical reasons this picture 
goes under the name of Replica Symmetry Breaking (RSB). 

Since analytic progress on finite dimensional system is hard, the study of 3D systems 
has mainly been based on numerical simulations. Most of the numerical simulations have 
investigated systems of size up to V = L 3 = 16 3 and temperature values greater or equal 
than half of the critical temperature (T c ): the results are in very good agreement with the 
RSB picture. Using the best techniques that are available to us today it is impossible to 
thermalize in a reasonable amount of time systems of this size at much lower temperature: 
because of that at today the very low temperature region cannot be explored using Monte 
Carlo simulations. On the other end a different class of algorithms exists which allow the 
determination of the ground state at zero temperature for systems of comparable size : 
it is clear that it would be very interesting to compare the results of the RSB approach with 
the numerical results at zero temperature. 

We face however a difficulty: in the RSB approach one computes quantities at finite 
temperature in the infinite volume limit, while if we look first for the ground state in a finite 
volume and only later we send the volume to infinity we are exchanging the order of the 
two limits. One can assume in a tentative way that the two limits can be exchanged, in the 
sense that in the zero temperature limit the probability distributions of the free energies of 
different equilibrium states becomes the probability distributions of the local minima of the 
energy: this has been the point of view of the authors that have employed this technique 
in the last months, it is not incompatible with any of the findings obtained to date, and we 
will do the same in the following. 

Of course at zero temperature the ground state (i.e. the global minimum of the energy) 
of a finite volume system with quenched random couplings selected under a continuous 
distribution is unique, so that interesting information can be extracted for example if we 
consider the effects of an external perturbation. 
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Different choices are possible. In a recent, very interesting paper, Palassini and Young 
[/J have shown that it is possible to get useful information about the nature of the low T 
phase of 3D Ising spin glasses with Gaussian distributed couplings by studying the behavior 
of the ground state after changing the boundary conditions of an L 3 lattice system from 
periodic (PBC) to anti-periodic (ABC). They analyze ground states (GS) obtained with 
the same realization of the quenched disorder and different boundary conditions (or with 
disorder realizations that differ only in some locations) for L < 10. In this paper we perform 
a more detailed analysis of the same systems by considering a larger set of observables. Here 
we use a modified genetic algorithm which is more efficient than the original algorithm and 
that will be described in details elsewhere ||. This modified algorithm has allowed us to 
study larger systems (L < 14) than in former work. We arrive at the conclusion that the 
whole set of data strongly suggests that the picture based on Replica Symmetry Breaking 
is correct in describing the behavior of 3D Spin Glasses. 



As we have already discussed in the introduction we consider a 3D Ising spin glass with 
quenched couplings assigned under a Gaussian probability distribution with zero expectation 
value, and we study the behavior of the ground state after changing the boundary conditions 
(BC) of a simple cubic lattice system of size V = L 3 from periodic (P) to anti-periodic (AP). 
When we do such a change the ground state is usually modified and in the new ground state 
some spins are reversed. Let us call a(i) and r(i) the spins in the ground states with periodic 
and anti-periodic boundary conditions respectively. The local overlap on site i is defined as 



II. SOME THEORETICAL CONSIDERATIONS 



A. Four Possibilities 




a(i) r(i) . 



(1) 



The link overlap is defined on the links and it is given by 



qi{i,H) = q(i) q{i + n) 



(2) 
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where /i is the direction of the link, that can take D positive and D negative values. If 
two spin configurations differ by a global reverse of the spins their link overlap is identically 
equal to one (while their overlap is equal to —1). We will consider in the following the case 
where we flip the boundary conditions in the direction x, while we leave them unchanged in 
directions y and z. 

In a first approximation (neglecting the points at the boundaries, see next subsection) 
we can define the interface among the flipped and the non flipped region as the sets of link 
where qi(i,fi) = —1. We are interested in finding out the geometrical properties of this 
interface in the infinite volume limit. Of course the probability of finding the interface on a 
random link is given by 

p=\{l-qi), (3) 

where qi is the disorder expectation value of g/(i,/i), averaged over sites % and directions jj,. 
Here we discuss a few possible situations: 

1. The interface is confined in a region of width L z (with z < 1): inside this region 
the interface may have overhangs. The interface density goes to zero as L~ a with 
a > 1 — z. This is what happens in ferromagnetic models, both with ordered and 
disordered Hamiltonians, where a — 1, al least in the ordered case. We will see 
immediately that this possibility is completely excluded by the data. 

2. The wandering exponent z is equal to one and the interface density goes to zero as 
L~ a . We also assume than in the large volume limit the interface is a fractal object 
with fractal dimension D s = D — a (where the space dimension D is equal to 3). More 
precisely we assume that if we define a continuous coordinate (in the interval [0 — 1]) 
as in the infinite volume limit the interface becomes a fractal object defined on the 
continuum characterized by a fractal dimension D s . In other words we are assuming 
that the interface is not a multi-fractal (i.e. that it can be characterized by a single 
fractal dimensions as usually happens for fractal objects one defines on a lattice in 
statistical mechanics), and so the relation D s = D — a is consistent. 
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3. The exponent a is equal to zero and the density p goes to a non zero value (i.e. qi 
does not go to 1 in the large volume limit). Here the interface is space filling and in 
the above described continuum limit, the interface is a dense set of measure 1. We 
expect that the probability that the interface does not intersect a region TZ, whose size 
is proportional to the system size, goes to zero when the volume goes to infinity. On 
the contrary in the previous case such a probability would be a non trivial function of 
1Z. As we shall see later, this possibility is the one realized in the case of RSB. 

4. There is a last possibility, which is unusual, and we mention for completeness (al- 
though its description will take a disproportionate amount of space): this is when the 
wandering exponent z is equal to one and a is greater than zero, but in continuum 
limit the interface does not become a fractal object of dimension D s = D — a, but it 
becomes a dense set. This happens for example if we consider a set of parallel planes 
at distance proportional to L a or a set of isolated points at distance Ln. In both cases 
the density goes to zero as L~ a , but the resulting object is not a fractal in the infinite 
volume limit. A similar case happens in three dimensions if we take a random walk of 
length £ = L 3 ~ a with a < 1. The system is a fractal of dimension D s up to a distance 
£(L), while it looks homogeneous at larger distances. 

In this case the density- density correlation function will scale as 



in the region 1 C i <C £{L), where C s = D — D s is the co-dimension of the fractal (the 
underlying random walk in our case). We will denote by the upper bar the average 
over the quenched disorder. On the contrary in the region £(L) <C x the correlation 
function will scale as 



equation (f|) implies that C s > a. Moreover the two formulas must be consistent in the 



C(x,L) = p(x)p(0) (xL~ a x 



-Cs 



(4) 




(5) 



The condition that at large distances the term in equation (|5|) dominates over the one in 
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crossover region where x £{L). Imposing this condition one finds that L Cs oc 

L~ 2a and therefore we obtain that 

t(L)<xL», U = £-<1. (6) 

This case is anomalous in the sense that usual scaling is not valid and there is a 
crossover length which increases as a fractional power of the side of the system. The 
usual scaling law (valid in the region x = O(L)) 

C(x,L) = L- 2a f(xL- 1 ) (7) 

is not satisfied and it is replaced by the unusual relation 

C(x, L) = L- a ~ c °f(xL- 1 ) + AL~ 2a g{xL- 1 ) . (8) 

Only in the case a = we are in a familiar situation: in this case the crossover length 
£ does not diverge for large L and we recover the third case of our list. 

B. The Effect of Changing Boundaries 

In the replica approach one finds that at zero temperature there are many different local 
minima of the Hamiltonian such that even in the limit where L goes to infinity the difference 
among the energy of the ground state and the energy of such minima is of order 1. It is 
crucial that the among these minima there are some whose site overlap q and link overlap 
qi with the ground state remains different from 1 in large volume limit. It is also important 
to assume that for large volumes there exists a function / such that qi = f(q 2 ), in other 
words that for fixed overlap q the link overlap qi does not fluctuate when the volume goes 
to infinity. 

When in a 3D spin glass the boundary conditions in the x direction are flipped from 
positive to negative it is known that the total energy changes by an amount AE, which 
increases as a power of L. In a ferromagnet, where the ground state is unique, the ground 
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state obtained under anti-periodic condition would be locally similar to the ground state 
obtained under periodic boundary conditions. In other words in any region of side M, if 
the region does not intersect the interface, which is just a flat surface, the ground state spin 
configuration obtained under anti-periodic boundary conditions will be equal (or equal to 
the reverse) to the ground state obtained under periodic boundary conditions. 

A similar conclusion holds in the case of a spin glass behaving according to the predictions 
of the droplet model [f|, where the boundary may be a corrugated surface: in this case we 
expect that the link overlap among the ground states obtained with PBC and ABC tends 
to one in the infinite volume limit. On the contrary in the case of a RSB like behavior the 
ground state obtained under ABC can be locally similar to one of the low energy excited 
states: in this case the link overlap among the two ground states can tend to a value different 
from one as L — > oo. 

That an excited state may be selected as new ground state when changing the boundary 
conditions is strongly suggested by the fact that the energy difference among ground states 
obtained with periodic and anti-periodic boundary conditions increases with L, while the 
difference among the ground state energy and the excited state energy at fixed boundary 
conditions remains fixed. 

The nature of the excited state that is selected when we change the boundary conditions 
cannot easily derived using simple arguments. It is reasonable to assume that in three 
dimensions, where AE which increases as a power of L, the new ground state will be as 
different as possible from the old ground state: one expects the overlap q among the two 
ground states to vanish in this limit. 

The interest in changing the boundary conditions is also due to the fact that this is also 
a convenient method to to study the properties of the excited states. Other methods can be 
used to clarify similar questions, but we will not use them here. 
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C. Gauge Invar iance 



In order to change the boundary conditions in the x direction we can change the sign 
of all the couplings connecting the plane x = x with the plane x — x + 1. The choice of 
xq does not matter: the ground state energy does not depend on x , and the ground state 
obtained after a different choice of the reversed plane (e.g. x = Xi) can be obtained starting 
from the ground state where anti-periodic boundary conditions have been enforced at xq, 
by flipping all the spins in the interval Xq < x < x\. 

In other words going from periodic to anti-periodic boundary conditions is a global 
change that locally is not visible. At this end it is convenient to define quantities which are 
gauge invariant, in the sense that they do not change when the plane where the anti-periodic 
boundary conditions are imposed. 

Let us consider a very simple example: the correct definition of qi(i, +x) is 

qi(i, +x) = q(i) q(i + x) for x ^ x , 

qi(i, +x) = —q(i) q{i + x) for x = x , (9) 

where x is the first component of the three dimensional vector i. 

A more complex case is the definition of the window overlap [|H]] in a box of side M 
which intersects the plane x = xq. Here the overlap qu can be defined as 

E (io) 

where A and B are the set of points of the box which are respectively on the right and on 
the left of the plane x = xq. 

The advantage of this prescription (which can be easily generalized to more complex 
situations) is that in absence of a magnetic field nothing depends on the actual position of 
the plane where the boundary conditions are enforced: the fact that the reversal of the spins 
is imposed at a given value of x is immaterial. The presence of this translational invariance 
is very effective from a practical point of view: the expectation value of qi does not depend 
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on the point, so that data can be taken on the whole lattice. The possibility of doing 
measurements on the whole lattice and not only far from the point where the boundary 
conditions are changed reduces the statistical errors in a substantial way. 

III. A FIRST LOOK TO THE NUMERICAL RESULTS 

We have computed 16000 couples of ground states on three dimensional lattices of linear 
size L = 4, 73740 for L = 6, 8532 for L = 8, 2434 for L = 10, 910 for L = 12 and 670 for 
L = 14: for each sample of the quenched disorder (Gaussian couplings) we have computed 
the two ground states under periodic and anti-periodic boundary conditions, using a modified 
genetic algorithm which will be described elsewhere || . 

A. The Average Link Overlap 

The introduction of anti-periodic boundary conditions breaks the discrete rotational 
invariance of the simple cubic lattice, so that the expectation value of qi(i,x), i.e. of the 
link overlap in the x direction, can be different from that of the link overlap in the other 
two directions. In the ferromagnetic case (without disorder) the change of the boundary 
conditions produces an interface in one of the y — z planes. We will call perpendicular 
link overlap qp the expectation value of qi(i,x) averaged over the lattice sites, since it is 
perpendicular to the y — z interface, while we will call transverse link overlap q T the average 
over sites of the quantity | (qi(i, y) + qi(i, z)). 

In the ferromagnetic case (without disorder) one finds that 

l- gp =I l-g T = 0. (11) 

We show in figure ([!]) the numerical data for 1 — q? and 1 — qp in our 3D spin glass. It is 
clear that the difference qr — qp goes to zero with r: it can be very well fitted as aL~ b with 
b ~ 2.5 and a normalized x 2 °f order one. 
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It is also evident that 1 — qp is with a very good approximation a linear function of L^ 1 
in the same way it would be in the ferromagnetic case, with the difference that here the 
extrapolation in the L —>■ oo limit is clearly different from zero. 

A closer look shows that 1 — qp is not exactly a linear function of L~ l , and it can be very 
well fitted (with a normalized x 2 °f order one) as a second order polynomial in L^ 1 . The 
same polynomial fit works very well also for 1 — qp. The two fits extrapolate to the same 
value (in the limit given by the statistical error). In the same picture we show also the data 
for 1 — qp — i.e. the value we measure for 1 — qp in our spin glass minus its a priori lower 
bound, which coincides with the value in the ferromagnetic case, and we note it has a very 
weak dependence over L. 



0.6 i- 
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0.2 - 



0. 1 - 1 - q transverse i — i — i 

1 - q perpendicular i — - * — ' 
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pure power law best fit to ( 1 - q transverse) 
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FIG. 1. 1 - q T , 1 - q P and 1 - q P - ~ versus ~ for L = 4, 6, 8, 10, 12, 14. We plot three 
polynomial best fits (of second degree in -^) to these three quantities and a simple power best fit 
to 1 — qp. 

A power law fit to a zero asymptotic value (i.e. to the form AL~ a ) of I — qp does not give 
acceptable results, while the data for l — qp can be fitted by a pure power law, with a ~ 0.27. 
A power law fits also works (even if with a higher value of x 2 ) f° r 1 ~ Qi(L) = 1 — |gr ~ \qp 
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with a similar value of a. 

It clear that it is not acceptable to fit the values of the two link overlaps with two 
different functional forms. Our conclusion is that the data do not support the droplet model 
prediction that q\ is zero in the large L asymptotic limit, and they hint for 

lim (1 - qi) = 0.245 ± 0.015 . (12) 

This estimate of the errors relays on the polynomial fit in powers of \. Considering a 
third order polynomial for qT does not change the situation. If we fit 1 — qr with a four 
parameter function of the form A + BL~ a + CL~ 2a we find an extrapolated value of A = qi = 
0.20 ± 0.06, which is consistent with the previous estimate: it is clear that in providing a 
complete estimate of the error over the extrapolated value one has to include the systematic 
uncertainty due to the possibility of using different functional forms for the extrapolating 
function. On the other end corrections proportional to a power of \ arise naturally, as it 
will be seen in the section on the correlations functions. 

The analysis of qx and qp shows that the RSB picture gives a consistent and satisfactory 
explanation of the numerical data. We will see in the rest of this note that this conclusion 
is supported by the analysis of many other quantities. 

B. The Hole Distribution 

It is interesting to consider a region of shape 1Z and to define the probability Pi (1Z, L) 
that a box of such a shape does not intersect the interface. In the case of a box of shape 
2 x 1 x 1 it is clear that PiiJZ, L) coincides with ~(1 + q{). 

Following Palassini and Young |7j we define 1 — P\(M,L) as the probability that the 
interface does hit a cubic box of side M, i.e. we take TZ = M x M x M. In figure (|2|) we 
show our numerical data for various values of L and M. It is possible to fit these data (if 
we restrict ourself in the region M < |-, in order to have comparable finite volume effects 
for the different L values) by a simple power of L (i.e. by the form AL~ a ): however the 
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power a of the best fits strongly depends on M (it is given by 0.32, 0.23 and 0.14 for M = 
2, 3 and 4 respectively), while in a meaningful fit we would expect to find a M-independent 
exponent. 

Second order polynomial fits in powers of -j- work well in the region M < ~ : they give 
for the extrapolated (1 — Pi (M, oo)) the values 0.35, 0.56 and 0.65 for M = 2, 3 and 4 
respectively. The data are consistent with the possibility that 1 — Pi(M, oo) goes to 1 when 
M —* oo. Indeed the extrapolated data reported before can be fitted as 1 — Pi(M, oo) = 
1 — AM' 1 with 7 close to one. 

It is clear that also in this case the data do not favor the kind of behavior predicted 
by the droplet model, where 1 — Pi(M, oo) = 0: they are far better consistent with the 
possibility that 1 — Pi(M, oo) is a non trivial function of M, which goes to 1 (and it is not 
zero) when M goes to infinity. 




FIG. 2. 1 - Pi(M,L) versus f for M = 2, min (L, 8), and L = 4, 6, 8, 10, 12, 14 



This conclusion is strengthened if we look for example to the plot of Pi (2, L) — Pi (3, L) 
and Pi(3,L) — Pi (4, L): the data for those two quantities are shown in figure (H). The 
droplet model suggestion that these quantities go to zero as a power of L when L goes to 
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infinity does not seem very convincing. 
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FIG. 3. The quantities Pi(2,L) - Pi(3,L) and Pi(3,L) - Pi(4,L) versus I for L = 4, 6, 8, 10, 



12, 14. 



IV. SCALING BEHAVIOR 

A. The Hole Distribution Again 

As we have discussed in the introduction if in the infinite volume limit the interface 
becomes a fractal, characterized by an unique fractal dimension D s , and the density of the 
interface goes to zero as L~ a , with a = D — D s , the probability of finding a hole in a box of 
size 1Z goes to a non trivial function of 1Z if we keep constant the size of the box in units of L. 
In the opposite case the interface would be become a space filling dense object. If we consider 
boxes of side M the previous argument shows that the probability that the box does intersect 
the interface 1 — Pi(M,L) should be a function of ^f. A glance to figure (||D shows that 
our numerical data do not exhibit such scaling: the interface is not a fractal characterized 
by a simple fractal dimension. On the contrary if we plot L~ 13 ln(Pi(M, L)) versus ^ (see 
figure (U) we find a reasonable scaling behavior, indicating again that 1 — P\(M,L) goes 
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exponentially to 1 when L goes to infinity at fixed ^f. 

The same conclusion holds very clearly if we consider the case of boxes of size L x L x 1 , 
i.e. planes, oriented in the y — z directions, i.e. parallel to the interface. The probability 
that such a plane does not intersect the interface P p (L) would tend to 1 in the situation 
(1) described in section ( [11 A| ), while it would tend to a non zero number in situation (2), 
The numerical data (in figure (|5[)) decrease very rapidly as a function of L (they can be 
reasonably fitted as aL~ s with 5 of order two, but a better fit has the form aexp(— bL c ): 
the asymptotic limit is zero in both cases). The first scenario we have presented in section 
( [HAD i.e. the one of an exponent z < 1, would imply that this probability goes to one 
asymptotically: this is certainly not the case. 

This set of data shows that the interface is space filling on a scale proportional to the side 
L of the system. Of course one could suggest that the numerical data are affected by very 
strong corrections to scaling: however it is not clear why this should happen (the scaling in 
figure (]4]) is very good). In order to see if we are able to detect any trace of these corrections 
it is convenient to consider other quantities. 
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FIG. 4. L- 1 - 3 ln(Pi (M,L)) versus f- for M = 2, min (L, 8) and L = 4, 6, 8, 10, 12, 14. 
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B. Other Scaling Laws 



If we go back to our boxes of side M we can define the window overlap JITJ as the value 
of the overlap q restricted to one of these boxes. We can define Pm,l(\q\) as the probability 
distribution of the absolute value of such window overlap: the probability is symmetric, so 
that we can consider only non negative values of q. The probability that the interface does 
not intersect the cube of side M is simply given by Pm,l(1)- 

In the case where the two ground states obtained under different boundary conditions 
are as different as possible the quantity 

Q 2 M,L = J dqq 2 P M , L {q) (13) 

should go zero at fixed j- at large M. A simple possible scaling behavior is that for large 
M and L 

Qm,l * L- s f(ML- 1 ) . (14) 

We show the numerical data for q\[ L L b versus ML~ X figure (|) (with 5 = .32). The data 
for L > 4 show a beautiful and very accurate scaling behavior. 
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FIG. 5. P P (L) versus L^ 1 for L = 4, 6, 8, 10, 12, 14. 
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The data in figures (§) and @ strongly suggest that asymptotic scaling laws are valid 
with good accuracy already for the lattice and block sizes that we have considered: sub- 
leading corrections are small, and they do not seem to be of exceptionally large size. More- 
over we find that for windows that occupy a finite fraction of the entire system in the infinite 
volume limit the window overlap distribution seems to become a delta function S(q). 
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FIG. 6. q 2 MtL L s versus ML' 1 for M = 2, min (L, 8) and L = 4, 6, 8, 10, 12, 14. 

Up to now all of our numerical evidence seems in perfect agreement with the RSB 
predictions, and the first two possibilities we had discussed in section ( [11 A| ) are excluded. 
In order to exclude the more exotic last possibility, and to present further evidence for the 
correctness of the RSB picture, it is convenient to study the correlation functions, both of 
qi and of q. This will be done in the next section. 



V. CORRELATION FUNCTIONS 
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A. The Link Overlap Correlation Functions 

We will present here a detailed analysis of the correlation functions: this is important 
since they carry a large amount of information that is crucial to distinguish among the 
different possible behaviors. 

As we have already discussed if we consider the case where there is an interface in the 
y—z plane we must distinguish among correlations in the transverse and in the perpendicular 
directions. Moreover in the case of vector like quantities like q\ we must distinguish among 
correlations in the direction of the link and correlations in the directions perpendicular to 
that of the link. For simplicity here we will consider only one correlation, i.e. the transverse 
link overlap correlation in the direction perpendicular to the link. More precisely the link 
overlap correlation function we consider is defined as 

G(d, L) = — y, z] y) q t (x, y,z + d;y) + qi(x, y, z\ z) q t (x, y + d,z;z) . (15) 

With this definition the correlation functions are periodic functions of d with period L, and 
they are symmetric functions of d, so that only the case < d < -| is interesting. 
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FIG. 7. The correlation functions G(d, L) versus w = f g + j^z^) with A = 1.5 for d > and L 
= 4, 6, 8, 10, 12, 14. The solid curves are smooth interpolations to the data point, and they turn 
out to be remarkably linear and the lines corresponding to different sizes are nearly parallel. 

We show in figure (|7|) the link overlap correlation functions at distance greater than 
for different L values. We plot them as a function of the variable 

Hs + r^) A - < l6) 

with A = 1.5. The reasons for this choice of the dependent variable will be clearer later: 
right now we just notice that the points look remarkably linear. 

If we want to concentrate our attention on the d dependence of the correlation function 
(neglecting a possible constant value at large distance) we can consider the quantity 

AG(d + -,L) = G(d,L) -G(d+1,L) . (17) 

Simple scaling implies that the function 

f(w, L) = AG(d + 1 L) (d- 1 + (L — d)- 1 ) 1 ^ , (18) 

where w = dL' 1 , should be independent from L for large L. The numerical data are shown 
in figure (§|). A good scaling is obtained with the choice of the value A = 1.5. 
We expect that the correlation function can be fitted at large L and d as 

G{d, L) = B(L) (d- 1 + (L - d)- 1 ) ~ X + A{L) 2 . (19) 

The form that we have taken is the simplest one which enforces the periodicity. Of course 
different forms are possible. In order to consider also the (i-independent part of the correla- 
tions functions it is convenient to fit the whole set of data shown in figure ([/]) using equation 



(19), where for simplicity we fix A = 1.5. The fits are good (they would also be good for 
similar values of A) so that the whole data can be reconstructed by the knowledge of the 
parameters A(L) and B(L), which are shown in figure (P|). 
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FIG. 8. The quantity f(w,L) = AG(d + \, + (L — d)- x ) 1+x with A = 1.5 as a function 

of d. Here L = 4, 6, 8, 10, 12, 14. 
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FIG. 9. The two fit parameters A(L) and B(L) as function of L for L = 4, 6, 8, 10, 12, 14. We 
show both a a power fit and a polynomial fit in 

While A(L) displays some dependence on L, there is no hint that the quantity B(L) 
should go to zero when L — > oo. Certainly both of them do not go to zero as L~ a with 
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a ~ .3 as expected in the case of possibility (4), and they can be consistently extrapolated 
to a non-zero value. If we would insist on power law fits we would find a value of a a = 0.15 
for A(L) and a ridiculous value of as = 0.03 for B(L). 

It is remarkable that A(L), that should tend to the link overlap in the L — >• oo limit, 
tends to 0.242, to be compared with the value 0.245 ± 0.015 (see equation (|T2"D) that we 
have derived from a completely different analysis: the picture that emerges from these data 
give consistent support to a RSB physical scenario. Also it is worth to comment more on 
the power fit of figure (||): even if they are compatible with the data it looks clear they 
are only minimally consistent with the physical behavior shown by the measured points. 
Typically these fits stay basically constant for all measured L values, and than they forecast 
a sharp descend at very high L values: the polynomial fits show on the contrary a consistent 
behavior in the measured L range and in the L — > oo limit. 

The second scenario we have discussed is in complete variance with the behavior of the 
correlation function data. The exponent A should be .3, not 1.5, as indicated by the data. 
This conclusion is in perfect agreement with the results of the hole probability distribution. 

We are now ready to analyze the possible correctness of the fourth scenario we have 
discussed (the interface is a dense set in the continuum limit). The results of figure @ show 
that we should have D s = D — A « 1.5. The fact that the interface would have a fractal 
dimension less than 2 is strange, but not impossible: the interface could look like the surface 
of T-lymphocytes, which is full of microvilli, i.e. quasi one dimensional objects. 

However in this case the correlation function (and therefore the correlation functions at 
fixed distance) should go to zero as L~ a : the correlation functions at distances of order L 
should go to zero as L~ 2a . The scaling observed in figure (||) should not be there. The two 
exponents a a and ag should be both equal to 0.3, which we have already seen it is not true. 

It is evident that also the exotic possibility (4) is not compatible with the observed form 
of the link overlap correlations (which are perfectly compatible with simple scaling laws 
without anomalies), and the only possibilities is given by the behavior implied by the RSB 
scenario. 

20 



B. The Overlap Correlation Functions 



Here we consider the overlap correlation function. In this case we can define two inter- 
esting correlation functions: the transverse correlation function 

C T (d, L) = ^— q(?, y, z) q(x, y,z + d) + q(x, y, z) q(x, y + d,z) , (20) 

^ x,y,z 

and the perpendicular correlation function 

c P {d, l) = — J2 ^0' y> z ) <i( x + d > y, z ) s{x,x + d) , (21) 

where S(x, x + d) = ±1 is factor which implements gauge invariance: it is equal to —1 if the 
line which connects the points x and x + d crosses the plane where boundary conditions are 
changed. 

The two correlations functions are respectively periodic and anti-periodic functions of D 
with period L (C P (L/2,L) = 0). 
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q-q perp 
q-q perp 
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FIG. 10. The overlap-overlap correlation functions divided times the factor (4 + 7—7) as a 



0.32 



function of j. The upper dots are for the perpendicular correlation function divided an additional 



factor (1 — tt). Lines are best fits (see the text). 
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We notice that 



q P {L) = C P {l,L), q T (L) = C T (l,L). (22) 

Consequently an analysis of the combined d and L dependence of the overlap-overlap corre- 
lation functions can give information on the origin of the L dependence of the link overlap, 
which coincide with the overlap-overlap correlation function at distance 1. 

According to the previous analysis of the window overlap we expect that the two corre- 
lations functions to scale as 

C{d, L) = f{dL- l ){d- 1 + (L - d)- l Y & . (23) 

The plot of figure (|10D shows a very good scaling, and two parameter polynomial fits work 
very well. Transversal and perpendicular correlation functions go to the same limit when 
x — > 0. Since the perpendicular correlation function is identically zero in d = §, we also plot 
the perpendicular C divided times (1 — ^r), removing in this way the main ^ dependence, 
which comes from a kinematical effect: it goes very smoothly to the correct limit when 
x -> 0. 

The extrapolated values are 0.732 ± 0.008 for the transverse correlation function and 
0.722 ± 0.005 for the perpendicular one. This is well consistent with the estimate quoted 
before (and obtained by measuring a completely different quantity) for the link overlap, 
qi = 0.755 ±0.015. 

As a consistency check we plot in figure (|TT| ) the correlation function Ct{L/2, L) as 
function of L. The data can be fitted very well by simple power fit C L~ s , with 5 ~ 0.32. 

It is clear that we have identified a physical mechanism which naturally generates an 
2 dependence in qi (correlation functions do usually feel the size of the system): simple 
scaling laws are valid and the ground state structure is non trivial. On the contrary a 
trivial structure (in the sense of the droplet model) would imply that the quantities f(s) 
extrapolate at 1 when s goes to zero (which is evidently not true) or they do not satisfy a 
simple scaling and there is a crossover region of length £(L) which goes to infinity with L. 
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Both these possibilities are not compatible with the data which strongly suggested the RSB 
scenario. 
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FIG. 11. The correlation C T {L/2,L) as function L = 4, 6, 8, 10, 12, 14. 

VI. CONCLUSIONS 

We have studied the overlap among the two ground states obtained under periodic and 
anti-periodic boundary conditions for 3D lattice of linear sizes L up to L = 14. We have 
analyzed many different observables: the average of the link overlap, the probability that the 
a cube of size M or a plane hits the interface, the window overlap probability distribution 
and the correlation functions of the link overlap and the overlap. 

The droplet model in the usual form assumes that the interface density (which is one 
minus the link overlap) goes to zero for large L as L~ a , and that the interface is a fractal 
with dimension D s = D — a = 3 — a. A small subset of our numerical data is consistent 
with the possibility the interface density could go to zero as L~ a , but other data cannot be 
fitted as a simple power and strongly suggest that the interface density remains finite in the 
limit L — > oo. The possibility that D s = 3 — a is completely excluded by the data. Some 
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of the data are compatible with a = 0.3 and D s = 1.5, but we have shown that this exotic 
possibility is not compatible with the rest of data. 

The RSB scenario is perfectly compatible with the whole set of data. The appropri- 
ate scaling laws for all the variables are satisfied with amazing accuracy. The observed L 
dependence of the interface density can be easily explained if we consider the scaling law 
appropriate for the overlap overlap correlation function. We also want to notice the coher- 
ence of these results and the recent detection of large scale excitations in spin glass ground 



states 11.1 



The behavior we find connects smoothly with the results obtained by simulations at 
finite temperature 0. For example we find that the overlap-overlap correlation decays as 
the distance to a power 5 = .3, which is not very far from the value obtained from simulation 



at finite temperature. Other results in this direction can be found in [|i~3 |. 

We think that we have conclusively shown that the ground state structure of 3D Ising 
spin glasses with Gaussian quenched random couplings is not trivial. 
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